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NOMENCLATURE 


A Qw/Tw 

a w sound speed based on the wall temperature 

B q heat transfer parameter 

C wall-law parameter 

C\ sublayer parameter 

c p specific heat 

c c i turbulence model constant for e 

c € 2 turbulence model constant for e 

Cp turbulence model constant 

c\ turbulence model constant 

C2 turbulence model constant 

D (A 2 + Ry / 2 

d\ constant defined in equation ( 18 ) 

o?2 constant defined in equation ( 19 ) 

0J3 constant defined in equation (20) 

k turbulence kinetic energy 

l index defined in equation (11) 

M free-stream Mach number 

Mr Ur / Q*w 

m index defined in equation (11) 

n index defined in equation (11) 

Pk production of kinetic energy 

Pvt turbulence Prandtl number 

q heat flux 

q w heat flux at the wall 


R 2 c p T w / Pvt 

Reg Reynolds number based on the momentum 
thickness 
T temperature 

T w temperature at the wall 

U velocity 

U c transformed velocity 

U+ U c /Ur, dimensionless transformed velocity 

U r frictional velocity 

WD shear stress 

y vertical distance from the wall 

U T y/v w , dimensionless normal coordinate 
e turbulence energy dissipation 

k Van Karman constant 

kj* Van Karman constant for temperature 

lit turbulence viscosity 

i/t kinematic turbulence viscosity 

u w kinematic turbulence viscosity at the wall 

p density 

p w density at the wall 

a*. turbulence Prandtl number for k 

a € turbulence Prandtl number for e 

Ofy turbulence Prandtl number for <j> 

r k/c, turbulence time scale 

r w wall shear stress 

<j> generalized turbulence length-scale variable 

uj e/k 



VI 




■ v r%$ 


iii 
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SUMMARY 


A critical assessment is made of the closure co- 
efficients used for turbulence length scale in existing 
models of the transport equation, with reference to the 
extension of these models to compressible flow. It 
is shown that to satisfy the compressible “law of the 
wall,” the model coefficients must actually be func- 
tions of density gradients. The magnitude of the errors 
that result from neglecting this dependence on den- 
sity varies with the variable used to specify the length 
scale. Among the models investigated, the k—u model 
yields the best performance, although it is not com- 
pletely free from errors associated with density terms. 
Models designed to reduce the density-gradient effect 
to an insignificant level are proposed. 

INTRODUCTION 

Calculations of high-Mach-number turbulent flows 
have become a major challenge in computational fluid 
dynamics in recent years. Traditionally, models devel- 
oped for incompressible flows have been extended to 
compressible flows with little or no modification when 
mass-weighted dependent variables are used. Our in- 
vestigation centers on the “dissipation-transport” equa- 
tion and its role in predicting the compressible law 
of the wall. This equation, or some other equation 
implying a length scale for the turbulence, is used in 
“two-equation” eddy-viscosity transport models and in 
full Reynolds-stress transport models. 

The Van Driest compressible law of the wall is 
derived from inner-layer similarity arguments that lead 
to the “mixing-length” formulas for velocity and tem- 
perature (Bradshaw, 1977). 

dU/dy = ( r w /p)V 2 /Ky (1) 

dT/dy = -q/[pcp(Wp)V2 KT2/] ( 2 ) 

Assuming r = t w , integration of (1) and (2) 
yields the temperature and velocity profiles 

T = C\T W q w / (cpT w ) — PrtU 2 /(2c p ) (3) 

U+ = UJU T = -In y+ + C (4) 


where Prt is the turbulent Prandtl number and is equal 
to k/k rp. The coefficient C\ is necessary here because, 
strictly, equations (1) and (2) do not apply in the vis- 
cous sublayer region. The quantities U T and y + are 
defined with respect to the physical properties at the 
wall as 

U T = (r w /p w )'/ 2 (5) 

and 

y+ = yU T /u w (6) 

The coefficients Cj and C are assumed to have 
the following forms, recommended from a data fit by 
Bradshaw (1977). 

C\ = 1 C = 5.2 + 95M 2 + 30.7 B q + 226 B 2 (7) 

where M r = U T /a w and B q = qw/(pu)C p U T T w ). 
The Van Driest transformed velocity, U c , is defined by 

U c = f U {-?-) x / 2 dU (8) 

Jo pw 

Substituting p/p w = T w /T and equations (3) and 
(7) into equation (8) yields 

U c = R'/ 2 [sm~ (9) 

where R = 2 c p T w /Prt, A = q w /T w and D = 
(A 2 + R) 1 / 2 . 

It should be noted that the Van Driest compress- 
ible law of the wall has been shown to accurately repre- 
sent experimental boundary-layer velocity profiles cov- 
ering a wide range of Mach numbers and Reynolds 
numbers (Fernholz and Finley, 1980). 

Figure 1 shows computed velocity profiles at 
Reg = 10,000, transformed using equation (9). The 
computations were obtained from the Navier-Stokes 
calculations of a Mach 5 flow over an insulated flat 
plate. TXvo versions of low-Reynolds-number k — e 
models were used: Chien (1982) and Launder and 
Sharma (1974). Also shown is a line representing the 
compressible law of the wall (eq. (4)), for k = 0.41. 
It can be seen that the calculated velocity profiles ex- 
hibit a steeper gradient than the expected value of 1 /k. 
This behavior was also observed in a Reynolds-stress- 
model calculation based on a low-Reynolds-number 


version of the Launder-Reece-Rodi model (Launder 
et al., 1975; Launder and Shima, 1989). To eliminate 
the uncertainties associated with the low-ReynoIds- 
number parts of the models, further calculations were 
performed that were based on the standard higb- 
Reynolds-number version of the k — e model coupled 
with compressible-wall-function treatments (Viegas and 
Rubesin, 1985; Huang, 1990).* Figure 2 shows that at 
M = 0.1, the predicted profile follows the Coles law of 
the wall, but at M = 5, the calculation displays an ear- 
lier rise of the profile in the log-law region. 2 It should 
be noted that the displacement of the profiles between 
the two calculations is the result of differences in the 
treatment of the viscous wall region and is not relevant 
to the present discussion. The excessive slope of the 
velocity profile has also been observed by other inves- 
tigators, both with the k — e model (Aupoix, 1990) and 
with a Reynolds-stress model (Viegas, J. R., private 
communication). 

The above results clearly indicate that standard 
turbulence models do not accurately predict boundary 
layers at high Mach numbers, and since the defect oc- 
curs with both the k— e and the Reynolds-stress models 
it appears that the trouble lies in the e-equation. 

This project is supported by NASA grant 
NCC 2-610-SI. We are grateful to Mr. J. G. 
Marvin and Dr. J. Shang (Wright R & D Center) for 
their encouragement, to Dr. J. R. Viegas and Mr. M. W. 
Rubesin for their review of the current manuscript, and 
to Mr. M. W. Rubesin for his suggestion of the form 
of (f> in equation (11). 

ASSESSMENT OF CLOSURE 
COEFFICIENTS 


In this section, attention is focused on the inner 
layer of a constant-pressure boundary layer, where con- 
vective transport is negligible and the eddy-viscosity 
hypothesis is formally applied. The velocity profile is 
governed by 


±, 

dy {M dy ) 


( 10 ) 


! Thc model coefficients are * 0.09, c«j = 1.44, c e2 = 
1.92, 0k = 1 and cr e = 1.17. The value of a e is altered from 1.3 
to 1.17 to give a better fit to the incompressible law of the wall. 

2 Experimental evidence has suggested that the correction func- 
tion for the wake region is nearly the same function of the empir- 
ically chosen Reynolds number as at low speeds. 


where pt is the turbulent viscosity and, for a k — e 
two-equation model, is defined as pt = c^pk^/e. The 
constant is fixed at 0.09 to satisfy —uv/k = 0.3 in 
local equilibrium. 

To discuss the problem in a general manner, a new 
variable, <j>, is defined in relation to the dissipation of 
the turbulence energy, e. 


<}> = p n k Tn f l (11) 


Equation (11) allows the construction of a variety of 
length-scale equations by the proper choice of n, m, 
and l. 

In the standard format of a high-Reynolds-number 
two-equation model, the transport equations for k and 
<f>, respectively, can be written as 


and 


d pt dk 

d , ut 9<t>^ _ d> <b 

= c ' pPk k - Cipt l 


( 12 ) 

(13) 


where a k and are turbulent Prandtl numbers for k 
and <j>, respectively. is the generation term 


„ du ,ou s 2 

Pk = = 


(14) 


where -W = rfp — r w /p. 

In equation (13), cj and C 2 are dimensionless co- 
efficients and are related to the corresponding coeffi- 
cients of the original e-equation, c el and c e 2 , according 
to 


Cl = lc t \ + m (15) 

C 2 = lc e 2 + m (16) 


In the standard k — e model, c e 2 = 1.92; in the newer 
two-equation models, such as the k — lj and k — r 
models, c € 2 is adjusted to 1.8 in order to provide a 
better fit of the decay law of homogeneous isotropic 
turbulence. A constant value of c e \ is generally chosen 
by computer optimization. 

To ensure that the plot of £/+ against In y + has 
a slope of 1 /k in the logarithmic region, a unique re- 
lation must exist between k and the other coefficients. 
The derivation of this relation is similar to that for the 
incompressible case, but with the following density ef- 
fects taken into account. First, in this region the usual 
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incompressible assumption of dk/dy = 0 is replaced 
by dpk/dy = 0, because -JTu/k = constant and t = t w 
= -puv. Second, the dissipation relation e = U^/ny is 
replaced by e = p) 3 / 2 / ny, because the assump- 
tion in this region is that e = = -Tiv(dU/dy). 

Finally, the density can no longer be factored out of the 
diffusion terms: d/dy(pt •••) ^ (1 / p)d / dy(yt ■ . .)• 

By substituting these assumptions into equation (13) 
and performing some mathematical manipulations, one 
can obtain 


- «i) 


= 1 + P 

where 


. ydp y 2 d 2 p ydp 2 

\8i + d2 7e? + 




(17) 


dj = n — m — 21 + 31 2 + 2ml — 2 nl + cj — (18) 


. 3 ;_,_ 

c <2 = n — m — -l + cj — 

2 


(19) 


d 3 = (n-m- h)(n - m - - ]-) - ^-c x — ( 20 ) 

2 2 2 2 cry. 


It can be seen that in order to extend the models for in- 
compressible flows to compressible flows without hav- 
ing to adjust the closure coefficients, the second term 
on the right-hand side of equation (17) must be neg- 
ligible compared to unity. This can be accomplished 
either with small values of d\, d 2 , and d 3 or with small 
density-gradient terms, ( yfp)dp/dy , —(y 2 /p)d 2 p/dy 2 , 
and (y 2 /p 2 ){dp/dy) 2 . 

Table 1 lists the coefficients d\ to d 3 for differ- 
ent turbulence models. The magnitude of the density- 
gradient terms associated with d\ to d 3 can be obtained 
by assuming that equations (3) and (4) are valid in the 
range = 30 to 1000 at the Reynolds number of 


Table 1. Coefficients associated with density-gradient 
terms 


Model 

n 

m 

1 

dl/l 2 

d 2 /l 2 

di/l l 

k — e 

0 

0 

1 

2.68 

0.18 

0.48 

k — 

0 

-2 

2 

0.63 

-0.13 

0.19 

k — cj 

0 

-1 

1 

0.56 

0.06 

-0.33 

a 

k — r 

0 

1 

-1 

2.47 

0.97 

0.70 


“The r-equation (r = k/e) (Speziale et ah, 1990) written in the 
form of equation (13) will lead to <j t = —1.44. 


the calculation in figure 1. The density-gradient terms 
are shown in figure 3 for M = 5 and an insulated 
wall. It can be seen clearly that the density effects at 
the Reynolds number of the computation are not neg- 
ligible and have to be taken into account in order to 
correctly predict the law-of-the-wall behavior. 

Equations (3), (10), (12), and (13) have been solved 
between y+ w 50 and y^~ ra 1000 for a Mach 5 flow 
over an insulated wall. The temperature profile is given 
by equation ( 3 ), with the free-stream temperature given 
at 15°C. The velocities at two boundaries are given, 
one fixed at the free-stream condition and the other cal- 
culated from the law of the wall (eq. (4)). The bound- 
ary conditions for k and 4> are estimated from k = 
(r w / P )/c]/ 2 and * = ( P ’>-xn-”'T« /2+ ”'/c'? /2 )/ 

(*y) 1 - 

Figure 4 shows the calculated density-gradient 
terms obtained using the k — e model. These results 
should be compared with those of figure 3, in which 
the profiles were obtained under the assumption that the 
compressible law of the wall has been satisfied. The 
square root of the right-hand side of equation (17) can 
be viewed as the ratio of the theoretical k (= 0.41) 
to the calculated k. This ratio is shown in figure 5. 
The value can be as high as 1.6 in the k — c model; in 
contrast, the k — u (u; = e/k) model of Wilcox (1988; 
also Wilcox and Traci, 1976; and Wilcox, D. C., pri- 
vate communication) displays only a mild increase of 
1 /k. This is because the coefficients d\ to d 3 associ- 
ated with density-gradient terms are all relatively small 
in the fc — tv model, as shown in table 1. Figures 6 (a) 
and 6 (b) show the calculated velocity profiles for the 
k — e and k — u models, respectively. These profiles 
agree with the behavior of k shown in figure 5. 

POSSIBLE REMEDIES 


One obvious way to predict the correct Van Driest 
law-of-the-wall profile is to incorporate the density- 
gradient terms directly into the coefficient c\? This 
is done by an iterative procedure and has been found 
to return the k — e model to the expected law-of-the- 
wall profile. As shown in figure 7, the value of c\ 
required in the e-equation drops to 50% of its original 
value at y + = 1000. This is too drastic a change to 
be acceptable. 

3 It appears (hat c t is a better candidate to adjust than a $ because 
equation (17) is derived under the assumption that cr $ is a constant. 
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Another feasible remedy is to find a new model 
that will make d\, d 2 , and d 3 all zero. Unfortunately, 
substituting equation (15) into equations (18) through 
(20) gives the solution n = m = l = 0. Since this is 
not a meaningful solution, it seems best to choose a 
model that will allow two leading coefficients to equal 
zero, such as d\ and d 2 — 0. Assuming n = 0 and 
m = — 1, the solution l = 5/6 and = 0.25 

makes both d\ and d 2 equal to zero. On the other hand, 
if one assumes n = 0 and m = —l, the requirement for 
<fj and d 2 to equal zero will lead to a k — \/u> model. 
These models give rise to an additional constraint on 
the relation between < 7 *. and < 74 , because c\ is a 

known value obtainable from the solution satisfying d\ 
and d 2 equal to zero. Calculations have demonstrated 
that these models do indeed reduce the density effects 
to insignificant levels: 0.43% and 0.15% for k — y/u 
and k — e 5 / 6 /fc models, respectively, at y + = 1000. 

CONCLUSIONS 

The present study has shown that the extension of 
incompressible turbulence models to compressible flow 
requires density corrections to the closure coefficients 
in order to satisfy the law of the wall (logarithmic law 
in Van Driest transformed coordinates). Equation (17) 
provides a way to estimate the error in boundary layer 
calculations. The k - uj model appears to be more at- 
tractive than the k-e model at high Mach numbers, be- 
cause the coefficients of the unwanted density-gradient 
terms are smaller. A length-scale transport equation 
can be devised to minimize the density effects and has 
proved successful, at least in boundary-layer flow. Fur- 
ther investigation of the proposed models in other flows 
is under way. 
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Mach = 5, Insulated Wall 



Figure 1. Transformed velocity profiles obtained by low-Reynolds number k — e models; Navier-Stokes 
calculations. 
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Insulated wall 



y + 

Figure 2. Transformed velocity profiles obtained by the high-Reynojds number k - e model with compressible 
wall functions; Navier-Stokes calculations. 
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M = 5, Insulated wall 



Figure 3. Density related terms calculated from the Van Driest compressible law of the wall. 
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